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Abstract. An outstanding issue in braneworld theory concerns the setting up of 
proper boundary conditions for the brane-bulk system. Boundary conditions (BC's) 
employing regulatory branes or demanding that the bulk metric be nonsingular have 
yet to be implemented in full generality. In this paper, we take a different route and 
specify boundary conditions directly on the brane thereby arriving at a local and closed 
system of equations (on the brane). We consider a one-parameter family of boundary 
conditions involving the anisotropic stress of the projection of the bulk Weyl tensor 
on the brane and derive an exact system of equations describing scalar cosmological 
perturbations on a generic braneworld with induced gravity. Depending upon our 
choice of boundary conditions, perturbations on the brane either grow moderately 
(region of stability) or rapidly (instability). In the instability region, the evolution of 
perturbations usually depends upon the scale: small scale perturbations grow much 
more rapidly than those on larger scales. This instability is caused by a peculiar 
gravitational interaction between dark radiation and matter on the brane. Generalizing 
the boundary conditions obtained by Koyama and Maartens, we find for the Dvali- 
Gabadadze-Porrati model an instability, which leads to a dramatic scale-dependence 
of the evolution of density perturbations in matter and dark radiation. A different 
set of BC's, however, leads to a more moderate and scale-independent growth of 
perturbations. For the mimicry braneworld, which expands like LCDM, this class 
of BC's can lead to an earlier epoch of structure formation. 
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1. Introduction 

A specific feature of braneworld models of gravity [U Ej and cosmology [31 HJ [5] is that 
the arena of observable events occurs in a lower- dimensional Lorentzian submanifold 
of the higher-dimensional space-time. For instance, in the scenario which we shall be 
considering, the brane is a four- dimensional boundary of a space-time with a noncompact 
('infinite') extra dimension, and an observer is in direct contact only with the induced 
metric on the brane. In this situation, one might worry about the well-posedness of 
physical problems from the observer's viewpoint. Indeed, phenomena taking place on 
the brane can depend upon physical conditions in the higher-dimensional bulk manifold 
which is inaccessible to an observer living on the brane. Consequently, the brane observer 
loses the ability of predicting classical physical phenomena on the brane. Indeed, it 
appears that no initial-value problem can be well posed on the brane since, no matter 
how detailed the specification of the initial conditions on a spacelike hypersurface on 
the brane, this does not uniquely determine events occurring either in the future or the 
past of this hypersurface. At best, one can uniquely predict the events in the future of 
a Cauchy hypersurface on the brane if one has the knowledge of the entire past to this 
hypersurface (and vice versa). This situation is illustrated by Fig. CD, which depicts the 
brane as a boundary of a bulk space with a noncompact extra dimension. Let a solution 
of the bulk-brane system be given. By smoothly perturbing data on a hypersurface S 
in the bulk at some distance from the brane, one modifies the solution on the brane only 
to the past of the hypersurface A and to the future of the hypersurface B\ hence, even 
complete knowledge of the data on the brane in the region between the hypersurfaces A 
and B is not sufficient for predicting the evolution to the past or future of this region. 
This is sometimes expressed by saying that the braneworld equations are nonlocal from 
the viewpoint of an observer living on the branejj] 

One might try to limit the space of solutions by specifying some appropriate 
boundary conditions for the full brane-bulk system. Much effort has been applied 
to formulate reasonable conditions in the bulk, by demanding that the bulk metric 
be nonsingular or by employing other (regulatory) branes. However, to date, these 
proposals have not been implemented in the braneworld theory in full generality - 
with certain success they were used only in a linearized theory under some additional 
approximations (see e.g., [SI CD, E]) — and it is clear that such an approach to boundary 

f Similar effects can also occur in our four-dimensional world. For an observer bound to the Earth, any 
local experiment has an uncertainty in its outcome because the events which lie outside the observer's 
past light cone are totally unknown but, in principle, can intervene the experiment at any moment 
of time. Thus, an unwanted cosmic gamma ray or gravity wave can suddenly affect the experimental 
situation. One can shield the experimental device from gamma rays, but not from gravity waves. 
Perhaps, we are fortunate that gravity waves (if they exist) are so feeble that even to detect them 
constitutes a big problem. In the opposite extreme case, were we permanently subject to gravity waves 
of strong amplitudes propagating in all directions, all our science of celestial motion would probably 
be useless (if not unattainable), just as the science of marine navigation becomes useless during a big 
storm. 
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Figure 1. Brane as a four-dimensional boundary of the five- 
dimensional bulk. 

By smoothly perturbing the data on a hypersurface 5 in the bulk at some 
distance from the brane, one modifies the solution on the brane only to the 
past of the hypersurface A and to the future of the hypersurface B; hence, even 
complete knowledge of the data on the brane in the spacetime region between 
the hypersurfaces A and B is not sufficient for predicting the evolution to the 
past or future of this region. 

conditions still leaves intact the property of nonlocality described above. 

The issue of cosmological perturbations on the brane in induced-gravity models 
remains poorly studied in the literature precisely because of the difficulties posed by 
the nonlocal character of this theory and the lack of proper boundary conditions. Thus, 
in treating the scalar perturbations on the brane, one is led to consider the problem in 
the bulk space using, for instance, the Mukohyama master variable [S] , which is a scalar 
function describing gravitational degrees of freedom in the bulk and obeying there a 
second-order partial differential equation. The corresponding initial-value problem was 
argued to be well posed in the brane-bulk system [ID] . However, in this case, one 
needs to specify initial conditions on the brane as well as in the bulk, with the obvious 
complication of having to deal with a function dependent on the coordinate of the extra 
dimension. This problem can be tackled to some extent in the linear approximation, 
but looks insurmountable when nonlinearities in the metric become significant — and, 
in principle, they can become significant in the bulk space much 'earlier' than on the 
brane. And — what is more important — one still has to face the problem of boundary 
conditions in the bulk, since not all possible values of the Mukohyama variable may lead 
to well-behaved solutions of the brane-bulk system. 

In an earlier paper [11], we suggested a new approach to the issue of boundary 
conditions for the brane-bulk system. From a broader perspective, boundary conditions 
can be regarded as any conditions which restrict the space of solutions. Our idea was 
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to specify such conditions directly on the brane which represents the observable world, 
so as to arrive at a local and closed system of equations on the brane. In fact, this is 
what some researches usually do in practice in the form of making various reasonable 
assumptions (see, e.g., [121 [13] f° r the issue of cosmological perturbations; and [H] for 
the case of spherically symmetric solutions on the brane). The behaviour of the metric 
in the bulk is of no further concern in this approach, since this metric is, for all practical 
purposes, unobservable directly. As first noted in [15], the nonlocality of the braneworld 
equations is connected with the dynamical properties of the bulk Weyl tensor projected 
on to the brane. It therefore seems logical to impose certain restrictions on this tensor 
in order to obtain a closed system of equations on the branejij] 

Neglecting the projected Weyl tensor on the brane [16J, though simple, is incorrect 
[7] since it is incompatible with the equation that follows from the Bianchi identity [see 
Eq. ([9]) below]. Hence, to proceed along this road, one needs to impose some other 
condition on this tensor. Perhaps, the simplest choice is to set to zero its (appropriately 
defined) anisotropic stress. This condition is fully compatible with all equations of the 
theory and eventually amounts to a brane universe described by a modified theory of 
gravity with an additional invisible component — the Weyl fluid, or dark radiation - 
having nontrivial dynamics. However, this is not a unique prescription and we shall 
discuss other possibilities in this paper applying some of them to the study of scalar 
cosmological perturbations on the brane. 

An important result of this paper is that our choice of boundary conditions (BC's) 
strongly influences the growth of density perturbations. For a class of BC's, a qualitative 
analysis of the linearized equations for scalar perturbations in a matter-dominated 
braneworld reveals the existence of an instability. In particular, the boundary condition 
which formally coincides with the approximate condition derived in [7] in the context 
of the Dvali-Gabadadze-Porrati (DGP) model [2], leads to a dramatic dependence 
of the growth of perturbations on spatial scale. In a matter-dominated braneworld, 
perturbations on small spatial scales grow many orders of magnitude larger than those 
on large spatial scales, indicating an early breakdown of linear regime. The cause of this 
instability is a peculiar gravitational interaction between matter and dark radiation on 
the brane which arises for this class of BC's. Our qualitative conclusions are supported 
by a numerical integration of the exact system of linearized equations. A different 
choice of BC's, however, leads to a perfectly 'normal' situation in which the growth of 
perturbations is moderate and virtually independent of scale. 

In the DGP model, the brane tension and the bulk cosmological constant are both 
set to zero, and the current acceleration of the universe expansion is explained as an 

| We admit that our current approach is quite radical because it effectively "freezes" certain degrees 
of freedom in the bulk; but its merit is that it apparently leads to a well-defined closed, local, causal, 
and, in principle, verifiable theory of gravity in four dimensions. If this proposal of specifying boundary 
conditions on the brane seems too radical, then one can view the results of this paper as an investigation 
of a certain class of approximations to the perturbation equations on the brane, some of which have 
appeared in the literature [3 [8] . 
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effect of extra dimension without severe fine tuning of the fundamental constants of the 
theory. However, an interesting situation arises also in the opposite case, where both 
these constants are large. As shown in a low-density braneworld exactly mimics the 
expansion properties of the LCDM model. In particular, a universe consisting solely of 
baryons with Q\> ~ 0.04 can mimic the LCDM cosmology with a much larger 'effective' 
value of the matter density Q m ~ 0.2—0.3. A conventional low-density universe runs 
into trouble with observations on account of its slow growth rate of perturbations. It is 
therefore interesting that such a model survives in the braneworld context which allows 
perturbations to grow much faster. 

Whether one can altogether do away with the notion of dark matter, as suggested 
in [T7|,[18] is a moot point, since, in addition to alleviating the so-called growth problem, 
dark matter also explains other observations including rotation curves of galaxies, etc. 
Whether all such observations can be accommodated in induced-gravity models is an 
interesting open subject which we set aside for future investigations. For recent results 
in frames of the the induced-gravity models, see our paper [T9] . 

This paper is organized as follows: In section [21 we introduce the theory and 
formulate the boundary conditions. In section [3j we apply the theory to the issue 
of scalar cosmological perturbations. We derive an exact convenient system of equation 
governing the growth of perturbations of pressureless matter and dark radiation and 
discuss boundary conditions of various type. In section HI we study numerically the 
evolution of perturbations in the DGP cosmology, and in section [5] we do this for the 
mimicry braneworld models. Finally, we discuss our results in section El 



2. Boundary conditions for the brane— bulk system 

We consider a generic braneworld model with the action given by the expression 
M 3 \f (K - 2A) - 2 / k] + f (m 2 R - 2cr) + f L (h ab , 0) . (1) 

. J bulk J brane J J brane ^ ' J brane 

Here, TZ is the scalar curvature of the metric g ab in the five-dimensional bulk, and R 
is the scalar curvature of the induced metric h ab = g ab — n a n b on the brane, where 
n a is the vector field of the inner unit normal to the brane, which is assumed to 
be a boundary of the bulk space, and the notation and conventions of [20] are used. 
The quantity K = h ab K ab is the trace of the symmetric tensor of extrinsic curvature 
K ab = h C aV c n b of the brane. The symbol L(h ab , <p) denotes the Lagrangian density of 
the four- dimensional matter fields whose dynamics is restricted to the brane so that 
they interact only with the induced metric h ab . All integrations over the bulk and brane 
are taken with the corresponding natural volume elements. The symbols M and m 
denote the five-dimensional and four-dimensional Planck masses, respectively, A is the 
bulk cosmological constant, and a is the brane tension. Note that the Hilbert-Einstein 
term f hrane m 2 R may be induced by quantum fluctuations of matter fields residing on 
the brane [21 H] . This idea was originally suggested by Sakharov in the context of his 
induced gravity model of general relativity [2T] . 
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The following famous cosmologies are related to important subclasses of action 

(i) The Randall-Sundrum model [1] is obtained after setting m = in JTJ. 

(ii) The Dvali-Gabadadze-Porrati model [2] corresponds to the special case where both 
the cosmological constant in the bulk and the brane tension vanish, i.e., A = and 
a = in flU). 

(iii) Finally, general relativity, leading to the LCDM cosmological model, is obtained 
after setting M = in (JTJ) . 

In addition to the above, braneworld cosmology described by (CEJ) is very rich in 
possibilities and leads to several new and interesting cosmological scenarios including 
the following [III 1221 123 I23| : 

• Braneworld universe can accelerate at late times, and the effective equation of 
state of dark energy can be to e g < — 1 as well as u> e g > — 1. The former provides 
an example of phantom cosmology without the latter's afflictions. 

• The current acceleration of the universe is a transient feature in a class of 
braneworld models. In such models, both the past and the future are described 
by matter-dominated expansion. 

• The braneworld scenario is flexible enough to allow a spatially flat universe to 
loiter in the past (z > few). During loitering, the universe expands more slowly 
than in LCDM (i.e., H(z)/H LCDM < 1), which leads to interesting observational 
possibilities for this scenario. 

• For large values of some of its parameters, the braneworld can mimic LCDM at low 
redshifts, so that H(z) ~ Hlcdm and w e s ~ —1 for z < few. At higher redshifts, 
the Hubble parameter departs from its value in concordance cosmology leading to 
important observational possibilities. 

It is therefore quite clear that the action (OQ) can give rise to cosmological 
models with quite definite attributes and properties. Confronting these models against 
observations becomes a challenging and meaningful exercise. At present, most tests of 
these models have been conducted under the assumption that the three-dimensional 
brane is homogeneous and isotropic [25] . Since braneworld cosmology has passed 
these tests successfully, the next step clearly is to probe it deeper by examining its 
perturbations. This shall form the focus of the present paper. 

The action (CQ) leads to the Einstein equation with cosmological constant in the 



bulk: 



g ab + Ag ab = , 



(2) 



while the field equation on the brane is 
m 2 G ab + ah ab = T ab + M 3 



(K ab - h ab K) , 



(3) 



Gravitational instability on the brane: the role of boundary conditions 



7 



where T a b is the stress-energy tensor on the brane stemming from the last term in action 
dl]). By using the Gauss-Codazzi identities and projecting the field equations onto the 
brane, one obtains the effective equation [TlJ [15] 

G a b + ^cfthab = 8nG c ffT a b + Q {-TTzQab ~ C a b) , (4) 



where 



(3 + 1 VM 6 



n 2am 2 

P = TT7^ ( 5 ) 



3M 6 

is a dimensionless parameter 




p + i 

is the effective cosmo logical constant, 

SttGW = JL-.-L (7) 
p + 1 m z 

is the effective gravitational constant, 

Qab = -EE a b — E ac E c b + - (^E cd E cd — -^E 2 ^j h a b (8) 

is a quadratic expression with respect to the 'bare' Einstein equation E a b = m 2 G a b — T ab 
on the brane, and E = h ab E a b- The symmetric traceless tensor C a b = n c n d C ac bd in (jlj) is 
a projection of the bulk Weyl tensor C a bcd- It is connected with the tensor Q a b through 
the equation following from the energy conservation and Bianchi identity: 



V a (Q ab 



-M 6 C ab )=0, (9) 



where V a denotes the covariant derivative on the brane associated with the induced 
metric h a b- 

The system of equations (j2J), ([3]) can have many formal solutions — too many to 
represent physically admissible cases. The issue that arises is what conditions should 
determine the space of its physical solutions. A common approach consists in demanding 
that the bulk space be singularity free. This approach is physically most appealing; 
however, because of obvious difficulties, it is very difficult to implement this idea in 
practice. (In the Appendix, we critically discuss the recent efforts [3 [8] to derive the 
proper boundary condition in this way.) Moreover, it may not be really necessary to 
impose this condition given that the entire observable world is represented by the brane 
only. Provided that the physical situation on the brane is regular, it may not be vital to 
worry about the situation in the bulk. This last consideration suggests that one might 
try to impose boundary conditions for the brane-bulk system not far away in the bulk 
but closer to home — directly on the brane or in a close neighbourhood of the brane - 
since this approach is likely to be the one most relevant for an observer residing in our 
(3+l)-dimensional universe. 

It is clear that the absence of the time derivatives of certain components of the 
traceless symmetric tensor C a b in Eq. results in a functional arbitrariness in the 
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dynamics of this tensor on the brane. Hence, it is reasonable to postulate boundary 
conditions in the form of certain constraints imposed on this tensor. In fact, conditions 
of this sort on the components of the tensor C ab have been applied in many papers 
as an approximation or a reasonable guess (see, e.g., [121 13] f° r a discussion in the 
Randall-Sundrum model). Here, we would like to consider such conditions as an exact 
physical principle. 

The tensor C ab can be decomposed through a convenient physically determined 
normalized timelike vector field u a on the brane (see |12|): 



where 



m 2 C ab = p C [u a U b + ^VabJ + 2W( a Ufe) + 7T ab , (10) 



q a b = h ab + u a u b (11) 

is the tensor of projection to the tangent subspace orthogonal to u a , and the covector 
v a and traceless symmetric tensor 7t ab are both orthogonal to u a . In this decomposition, 
the quantities pc, v a , and ir ab have the meaning of dark-radiation density, momentum 
transfer, and anisotropic stress, respectively, as measured by observers following the 
world lines tangent to u a . 

The anisotropic stress can be further conveniently decomposed into scalar, vector, 
and tensor parts in the following way: 

TTafc = D(aDb) ~ ^QabD 2 7l C + D {a 7l b) + TC ab . (12) 

Here, we have introduced the derivative operator D a acting on tensor fields T^jz tangent 
to the brane and orthogonal to u a according to the rule 

i)„r;t \\v„r;:!:::. (13) 

where the symbol T7 denotes the projection using q ab with respect to all tensorial indices. 
For example, 

D a Ti b = q c a q d b V c TT d . (14) 

If the vector field u a is hypersurface orthogonal, then q ab represents the induced metric in 
the family of hypersurfaces orthogonal to u a , and D a is the (unique) derivative operator 
on these hypersurfaces compatible with this induced metric. Furthermore, the covector 
field ir a and the symmetric tensor 7f ab in f ll2p are both orthogonal to u a and, for their 
unique specification, the latter can also be assumed to have zero trace and 'spatial 
divergence' (i.e., to be trace-free and 'transverse'): 

h ab Jr ab = , D a n ab = . (15) 

The first of these conditions implies D a n a = 0. 

There are no evolution equations for the tensor fields ir ab on the brane, which is 
a manifestation of the nonlocality problem discussed in the introduction. Therefore, 
boundary conditions on the brane can be specified by imposing additional conditions on 
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this tensor. The simplest possible choice of such conditions would be to set it entirely 
to zero: 



In this case, equation (Q gives the evolution equations for the remaining dark-radiation 
components pc and v a . However, other conditions on the components n c , n a , and 7f ab 
of the anisotropic stress are possible. Without introducing any additional dimensional 
parameters into the theory, one can relate the scalar 'spatial Laplacian' D 2 ttq to the 
energy density pc 



while setting n a = and 7f a b = 0. Here, A is some dimensionless covariant scalar. In 
this paper, we mostly assume A to be just a constant so that the boundary condition 
(j!7p forms a one-parameter family. 

Boundary conditions of the type ( 1T61) or ( TTTj) clearly depend on the choice of the 
timelike vector field u a used in the decompositions (flOj) and ([12!) . In order for these 
BC's to be covariantly specified, the vector field u a should be determined by the intrinsic 
geometric properties of the brane. One can think of several options as to the choice of this 
field. For example, one can identify u a with a timelike eigenvector (thereby demanding 
the existence of such an eigenvector) of a geometric quantity such as C a b, R a b, K a b, 
etc. We will not attempt to make a specific prescription in this paper. In studying 
cosmological scalar perturbations, it will be sufficient to assume that the unperturbed 
field u a is tangent to the world lines of isotropic observers — the only distinguished 
choice — and that it is perturbed linearly. In this case, condition ffl7|) represents a 
generalization of the approximate boundary condition derived in [7J E] for the DGP 
cosmology [2j [5] (there, the value A = —1/2 was obtained). 

The specific form ( Fl6l) or (TlTl) of our set of boundary conditions, of course, contains 
some arbitrariness. To us it seems to be the simplest one of the kind of boundary 
conditions that we propose (specified directly on the brane); however, we cannot justify 
its uniqueness at this moment. Therefore, it should be regarded as a proposal which 
demonstrates how the braneworld theory with our type of boundary conditions can 
operate in principle, and which can be tested observationally. Perhaps, one should 
keep in mind that the specific form of boundary conditions in the multi-dimensional 
space is part of any braneworld theory with large extra dimension, to be tested against 
observations. 

After the boundary conditions have been specified, the tensor C a b plays the role of 
the stress-energy tensor of an ideal fluid with equation of state like that of radiation 
but with a nontrivial dynamics described by Eq. (Q. In the literature, this tensor has 
been called 'Weyl fluid' [26] and, in the cosmological context, 'dark radiation' [27]. The 
stress-energy of this ideal fluid is not conserved due to the presence of the source term 
Qab m ©• Equation (0J together with boundary conditions such as those described by 
(JT0l)-( rr71) form a complete set of equations on the brane. 



TTafe = . 



(16) 



D\ c = Ap c , 



(17) 
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3. Scalar cosmological perturbations 



3.1. Main equations 

The unperturbed metric on the brane is described by the Robertson- Walker line element 
and brane expansion is described by pfl EJ [22] 



H 2 + - 2 

a 1 



P 



a 2 
3m 2 " + P 



P 



a 



3m 2 



A 
"6 



C 



where p is the matter density, 
„ 2m 2 



M 3 



;is) 



(19) 



defines a new fundamental length scale, k = 0, ±1 describe different possibilities for the 
spatial geometry and the term C/a A (C is a constant) is the homogeneous contribution 
from dark radiation. The two signs in (1181) describe two different branches corresponding 
to the two different ways in which a brane can be embedded in the Schwarzschild-anti- 
de Sitter bulk [HE]. In [22], we classified models with lower (upper) sign as Brane 1 
(Brane 2). Models with the upper sign can also be called self- accelerating because they 
lead to late-time cosmic acceleration even in the case of zero brane tension and bulk 
cosmological constant [5]. Throughout this paper, we consider the spatially flat case 
(k — 0) for simplicity. 

On a spatially homogeneous and isotropic brane, the timelike vector field u a used 
in the decomposition ({TO]) coincides with the four-velocity field of isotropic observers 
comoving with matter, and the homogeneous dark-radiation density pc is related to the 
constant C by p c = — ?>m 2 C/a 4 . Due to spatial homogeneity and isotropy, we have 
7i a = and 7i a b = 0, and equation ( fTTl) then implies pc = C = (for A ^ 0). 

Scalar metric perturbations are most conveniently described by the relativistic 
potentials $ and \l/ in the so-called longitudinal gauge: 

ds 2 = -(1 + 2<$>)dt 2 + a 2 (l - 2^/)dx 2 . (20) 

We denote the components of the linearly perturbed stress-energy tensor of matter 
in the coordinate basis as follows: 

/ ^(n + Sn\ -d iV \ 

(21) 

(p + Sp)^ 



rpO> 



13 



V 



<P + Sp) 

n2 



J 



where dp, 8p, and v are small quantities. Similarly, we introduce the scalar perturbations 
&Pc, v Ci an d 5n c of the tensor C a b in the coordinate basis: 

/ 3m 2 C ^ „ \ 



m 2 C a p 



- Sp c 

diV C 



-diV C 



m 2 C s 



5n l 



(22) 
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where, according to (fT21) . Sir^ = didjSnc — j$o~ijd 2 5iic, and the spatial indices are raised 
and lowered with 8 % y We call v and vc the momentum potentials for matter and dark 
radiation, respectively. 

Then, equation (j3J) together with the stress-energy conservation equation for matter 
and conservation equation ([9]) for dark radiation result in the following complete system 
of equations describing the evolution of scalar perturbations on the brane: 



* + 3(1 + 7)#i + [2H + 3H 2 (1 + 7) $ 

5p_ 

2m 2 



^V 2 * + t^V 2 ($ - *) 



a 



7 



2 

^A 



+ 1 + 



t5S 



l 2 X) 2m 2 ' 



1 



dp + 3H(Sp + <5p) = 3(p + p)V + -^V 2 v 



v + 3Hv = 5p + (p + p)$ , 

1 „ 12m 2 C • 
a 



5 p c + AH5p c = f^VV 



Vc + 3Hv c 

— V 2 * = 

a 1 

m 2 £ 2 X + H$ 



l r 4m 2 C" m 2 £ 2 A 

-<$Pc : — $ 

3 H a 4 3# 



Am 

2m 2 



m 2 £ 2 A 
3a 2 



V 2 ($ - 



Ay, 



£ 2 A/ 2m 2 



+ 



' ex 
1 + - 



m 2 £ 2 X ' 



5 Tlr 



m 2 £ 2 X 
4~~ 



l + 3 7 ) ($-*). 



(23) 

(24) 
(25) 

(26) 
(27) 
(28) 
(29) 

(30) 



Here, we use the following notation: S is the entropy density of the matter content 
of the universe, r = (dp/dS) , c 2 = {dp/ dp) s is the adiabatic sound velocity, the 



time-dependent functions A and 7 are given by 

p + a 



X = H Z 



7 



1 

3 



3m 2 
' HX 



a 



3m? 



A 
"6 



C 



p + p AC 



77T 



p + a 1 
1 — 

3m 2 £ 2 



A 
"6 



C 
771 



(31) 



(32) 



and the perturbations A m and Ac are defined as 

A m = 5p + 3Hv , A c = 5p c + 3Hv c . (33) 

The overdot, as usual, denotes the partial derivative with respect to the time t. 

The system of equations (|23|) - (|30|) generalizes the work of Deffayet [10] (for the 
DGP brane) to the case of a generic braneworld scenario described by (pQ), which allows 
non-zero values for the brane tension and bulk cosmo logical constant. It describes two 
dynamically coupled fluids: matter and dark radiation. It is important to emphasize 
that the evolution equations (|26|) . (|27|) for the dark-radiation component are not quite 
the same as those for ordinary radiation. Of special importance are the last two 
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source terms on the right-hand side of (1271) which lead to nonconservation of the dark- 
radiation density. Thus, the behaviour of this component is rather nontrivial, as will be 
demonstrated in the next section. 

It is also interesting to note that the perturbations in dark radiation formally 
decouple from those in ordinary matter in the important limiting case M — > 
(equivalently, £ — ► oo), for which the system fT23l) - fl30l) reproduces the corresponding 
equations of general relativity (after setting 7 = c 2 s ). 

From equations (|24l-(l29|). one can derive the following useful system for 
perturbations in pressureless matter and dark radiation in the important case C = : 

+ <*> 
v c + AHv c = 7 A C + ( 7 - I) A m + ^^^c , (35) 

Sp c + 4H5 Pc = \v 2 v c , (36) 



a 2 



where 



A = ^ (37) 
P 

is the conventional dimensionless variable describing matter perturbations. 
3.2. Boundary conditions 

The system of equations (T23T) — (T30]) . or fl34l) - fl36l) . describing scalar cosmological 
perturbations, is not closed on the brane since the quantity 5n c in (13"U|) or (1331) (hence, 
the difference $ — 1 J r ) is undetermined and, in principle, can be set arbitrarily from the 
brane viewpoint. This is a particular case of the nonclosure of the basic equations of 
the braneworld theory which we noted earlier. As discussed in the introduction, our 
approach to this issue will be different from that in other papers in that we will not 
look to the bulk to specify boundary conditions but, rather, will specify such conditions 
directly on the brane. 

A general family of boundary conditions on the brane is obtained by relating the 
quantities ttc and pc as discussed in the introduction. Equation (|T7|) implies C = in 
the case A 7^ and, in the linearized form, becomes 

\v 2 67r c = A5p c . (38) 

In most of this paper, A shall be assumed to be a constant. By virtue of fl30|) . this 
relates the difference $ — ^ between the gravitational potentials to the perturbation of 
the dark-radiation density Spc : 

1 A A 

-V 2 ($ - *) = -—-—^8 pc ■ (39) 

a 1 m l l l \{\ + 37) 

For the boundary condition (|38|) . one can derive a second-order differential equation 

for 5pc by substituting for vq from (1361) into (1331) : 

5p c + (10 - i)H5p c + 4 \H + 3(2 - -f)H 2 } 8p c 
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lv 2 

a 2 



l5pc + wh^ 5pc + ( 7 - 1) Ar 



(40) 



Equations (IM]) and (1401) will then form a closed system of two coupled second-order 
differential equations for A and 5pc- From the form of the right-hand side of ( l40l) one 
expects this system to have regions of stability as well as instability. Specifically, a 
necessary condition for stability on small spatial scales is that the sign of the coefficient 
of V 2 5pc on the right-hand side of ( l4TJi) be positive. This leads to the conditional 

A > -|-y(l + 37) ■ (41) 

From (I3"2j) , we find that 7^— 1/6 in a matter-dominated universe, and condition (|4"T1) 
simplifies to 

A>± (42) 

We consider two important subclasses of fl38|) which we call the minimal boundary 
condition and the Koyama-Maartens boundary condition, respectively: 

3.2.1. Minimal boundary condition. Our simplest condition ( jl~6l) corresponds to setting 
A = 0. Then, from (1301) we obtain the relation $ = the same as in general relativity. 
Under this condition, equations (I23l) -( l25l) constitute a complete system of equations 
for scalar cosmological perturbations on the brane in which initial conditions for the 
relativistic potential $, $ and matter perturbations 5p, 5p, v can be specified quite 
independently. Once a solution of this system is given, one can calculate all components 
of dark-radiation perturbations using (128]) and (1291) . Thus, with this boundary condition, 
equations (l26l) -( T29l can be regarded as auxiliary and can be used to felicitate and 
elucidate the dynamics described by the main system f[23]) — (|25]) . We should stress that 
only the quantities pertaining to the induced metric on the brane (<3>, \1/) and those 
pertaining to matter (5p, 5p, v) can be regarded as directly observable, while those 
describing dark radiation (Spc, t>c) are not directly observable. 

3.2.2. Koyama-Maartens boundary condition. In an important paper [7j, Koyama and 
Maartens arrived at condition fl38l) with 



A = -\- (43) 

This boundary condition was derived in [7] as an approximate relation in the DGP 
model valid only on small (subhorizon) spatial scales under the assumption of quasi- 
static behaviour. It was later re-derived in [5] under a similar approximation. However, 
according to the approach to boundary conditions taken in the present paper, one can 
regard (1431) as being a linearized version of the more general relation (|T7|) which is valid 

§ In the absence of matter on the brane, A m = 0, equation becomes a closed wave-like equation 
for the scalar mode of gravity, and condition (j¥T|) becomes the boundary of its stability domain. The 
existence of such a scalar gravitational mode is due to the presence of an extra dimension. 



Gravitational instability on the brane: the role of boundary conditions 



14 



on all spatial scales. We call it the Koyama-Maartens boundary condition, although one 
should be aware of the different status of this relation in our paper, where it is regarded 
as an exact additional relation, and in [7J [8] , where it is derived as an approximation to 
the more complicated situation. 

In discussing the small-scale approximation in quasi-static regime, it was argued in 
[7J that equation fl36|) permits one to neglect the perturbation t> c in fl35l) . which, together 
with (1381) . will then transform (|34"|) into a closed equation for matter perturbations!]]] 

A + 2HA = 6km & , 6km = 1 + - lMl±jl±i . (44) 



2m 2 ' A £2\ \2A+ §7(1 + 37) 

Some cosmological consequences of this approach are discussed in [28] . As in general 
relativity, this equation does not contain spatial derivatives; hence, the evolution of 
A is independent of the spatial scale. We would obtain equation of the type (jUJ) for 
perturbations had we followed the route of [7] or [8] in finding approximate solutions 
of perturbation equations in the bulk and using the quasi-static approximation (some 
difficulties with this approach are discussed in the appendix). However, in the approach 
adopted in this paper, as will be shown in the following two sections, numerical 
integration of the exact linearized system (!3~4"]) - (!36]) does not support approximation 
(1441) . No matter what initial conditions for dark radiation are set initially, one observes 
a strong dependence of the evolution of matter perturbations on the wave number. In 
particular, it is incorrect to neglect the quantity vc on small spatial scales, since it is 
precisely this quantity which is responsible for the dramatic growth of perturbations 
both in A m and in Spc on such scales. 

From equations (jUJ) and ( 1421) . we can see that the minimal and Koyama-Maartens 
boundary conditions generally lead to unstable evolution. This will be confirmed by 
numerical simulations in the next section. 

3.2.3. Scale-free boundary conditions. Evolution of perturbations in the stability region 
(|4T|) & (f42|) shows little dependence on spatial scale. It is interesting that there also 
exists an important class of boundary conditions leading to exact scale-independence. 
We call these, for simplicity, scale-free boundary conditions. To remove the dependence 
on wave number altogether and thereby obtain a theory in which perturbations in matter 
qualitatively evolve as in standard (post-recombination) cosmology, it suffices to set the 
right-hand side of (1351) identically zero: 

i ^ 2A i + 3 7 r 

which, in view of equation (1281) . can also be expressed in a form containing only the 
geometrical quantities Snc, Ac, and In this case, the perturbations 5vq and Spc 

|| Equation (I44|) was derived in [7] only for the DGP model and for the case A = —1/2; however, the 
argument can be extended to a general braneworld model and a general value of A in (|38|) . 



2 V^ttc = (1 - 3 7 )A m - 3 7 A C , (45) 
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in dark radiation decay very rapidly, according to equations (1551) and (|3"ET) . and ([511) 
reduces to the simple equation 

A + 2HA = 6 , e=(l + |2), (46) 

valid on all spatial scales. Equations ( |28l) and ( |30l) then lead to simple relations between 
the gravitational potentials $ and \I/ and matter perturbations: 

l V 2 $ = e— , -V 2 ^ = + — ) — . (47) 
a 2 2m 2 ' a 2 V £ 2 A7 2m 2 1 ; 

The difference $ — \1/ can be conveniently determined from 

> 2 (*-*) = ^„. (48) 

As can easily be seen from ( 1301) or (148p . the general- relativistic relation $ = ^ is 
not usually valid in braneworld models. An important exception to this rule is provided 
by the mimicry models discussed in section 15.31 

One can propose other conditions of type (14"5T) that lead to scale-independent 
behaviour. For instance, one can equate to zero the right-hand side of (T40|) . But it 
remains unclear how these may be generalized to the fully nonlinear case as was done 
in the case of the boundary conditions (138|) via equation (fTTj) . Nevertheless, in view of 
the interesting properties of scale-independence and the fact that perturbations in the 
stability region (j4Tl) & (j4*2l) behave in this manner, the consequences of ( 1461) need to be 
further explored, and we shall return to this important issue later on in this paper. 

Having described the system of linearized equations governing the evolution of 
scalar perturbations in pressureless matter and dark radiation, we now proceed to 
apply them to two important braneworld models: the popular DGP model (2j [5] and 
the 'mimicry' model suggested in [11]. It should be noted that these two models are 
complementary in the sense that the mimicry model arises for large values of the bulk 
cosmological constant A and brane tension a, whereas the DGP cosmology corresponds 
to the opposite situation A = and a = 0. 

4. Scalar perturbations in the DGP model 

Amongst alternatives to LCDM, the Dvali-Gabadadze-Porrati (DGP) model [2] stands 
out because of its stark simplicity. Like the cosmological constant which features 
in LCDM, the DGP model too has an extra parameter i = 2m 2 /M 3 , the length 
scale beyond which gravity effectively becomes five-dimensional. However, unlike 
the cosmological constant whose value must be extremely small in order to satisfy 
observations, the value i ~ cHq 1 , required to explain cosmic acceleration, can be 
obtained by a 'reasonable' value of the five-dimensional Planck mass M ~ 10 MeV. 
As pointed out earlier, DGP cosmology belongs to the class of induced gravity models 
which we examine and is obtained from (CD) after setting to zero the brane tension and 
the cosmological constant in the bulk (i.e., a = and A = 0). Under the additional 
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assumption of spatial flatness (k = 0) and C — 0, the modified Friedmann equation (fISI) 
becomes j5] 



2 2H p 



H z - 



3m 2 



(49) 



In a spatially flat universe, given the current value of the matter density and Hubble 
constant, £ ceases to be a free parameter and becomes related to the matter density by 
the following relation 

1 /i-a^ 2 



ftp 



t 2 Hl 



(50) 



which may be contrasted with = 1 — fl m in LCDM. 
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Figure 2. The DGP brane with the Koyama Maartens boundary 
condition A = —1/2. Growth of perturbations in matter (left) and dark 
radiation (right) on the DGP brane are shown for different values of the 
comoving wave number k/aoHo (indicated by numbers above the corresponding 
curves). The current matter density is chosen to be O m = 0.22, and the 
initial value of vq is set to zero. (Our results remain qualitatively the same for 
other values of the density parameter.) Note the dramatic ^-dependence in the 
growth of perturbations both in matter and in dark radiation. For comparison, 
A/a = 1 at these redsfnfts in the standard CDM model with Q m = 1. 



Linear perturbation equations for this model were discussed in [SJ El El ED]- An 
approximate boundary condition for scalar perturbations was obtained by Koyama and 
Maartens [7] on subhorizon scales, and it is described by equations (13"5j) . fl4"3"l) . For 



Gravitational instability on the brane: the role of boundary conditions 17 



convenience, we present system (!53|) - (j5o r |) for this casejj 

v c + AHv c = 7 A C + ( 7 - i) A m - ^^Spc , (52) 

5 pc + AH5 Pc = \v 2 v c ■ (53) 

In the DGP model, the general expressions (|3~T!) and (1321) for A and 7 in the case of 
pressureless matter reduce to 

x 2 / 1\ 11 

x = j{ H -j)' -' = 2im^W'e- (54) 




Figure 3. Shown is the logarithm of the ratio ()55f) for scalar perturbations in 
the DGP model described by the system of equations (f5TT) -(p)5 j) with comoving 
wave number h/a^Ho = 50. Initial conditions (at z = 10 3 ) are chosen such 
that vc = and vc = 0, so that the initial value of the ratio (f55j) is equal to 
zero. Observe that this ratio grows quickly and exceeds unity by several orders 
of magnitude at late times (low redshifts) demonstrating the breakdown of the 
quasi-static approximation. 

The results of a typical integration of the exact system of equations fl5TT) - fl53l) 
for different values of the wave number k are shown in figure [2J We observe a 
dramatic escalation in the growth of perturbations at moderate redshifts and a strong k- 
dependence for perturbations in matter as well as in dark-radiation (the y-axis is plotted 
in logarithmic units). These results do not support the approximation made in [7j, which 
assumes the left-hand side of equation fl52l) to be much smaller than individual terms 

If To facilitate comparison, we relate the notation of Ref. [7] with that of our paper. Our length scale 
I is related with the length scale r c of [7] by i — 2r c . Our quantity Sp is the same as in [7J, and our 
quantities v, Spc, vc, and Sttq in the DGP model are related to the similar quantities Sq, 5pE, &<1e, and 
Site of [7J by the equalities v = —adq, 5pc = —SpE, vq = aSqE, and Sire = — ci^Stte- 
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Figure 4. The DGP brane with the boundary condition A = 1/2. 

The parameters of the model are the same as in figure [2] but the y-axis is no 
longer plotted in logarithmic units. The two solid curves show the evolution of 
scalar perturbations corresponding to the comoving wave numbers fe/ao-ffo = 
(thin curve) and 40 (thick curve). Note that these two curves are almost 
indistinguishable which illustrates that the growth of perturbations is virtually 
scale-independent in this case. The dashed line shows the behaviour of scalar 
perturbations in the LCDM model. (In all cases O m = 0.22 is assumed.) 



on its right-hand side for sufficiently large values of k, and which leads, subsequently, to 
the scale-independent equation (JHJ. To show explicitly that this 'quasi-static' regime is 
rather unlikely, we have integrated the system (j5Tl) - (j5"51) for a sufficiently high comoving 
value k/aoHo = 50 and with the initial condition vc = and vc = 0, which sets both 
sides of equation ([52]) initially to zero. In figure EH we plot the logarithm of the ratio 



which is assumed to be < 1 in the quasi-static approximation of [7J. From figure [3] 
we see that this ratio soon becomes much larger than unity. This demonstrates that, 
even if vq and vc were small initially, during the course of expansion both quantities 
grow to large values casting doubt on the validity of the quasi-static approximation for 
small-wavelength modes. 

We would like to stress that our conclusions themselves are not based on the small- 
scale or quasi-static approximation. Indeed, we integrate the exact system of equations 
fl34l)-(l36l) on the brane, and the only ansatz that we set in this system is the boundary 
condition fl38l) . (|43l) . Our failure to find the regime in which the ratio (155]) remains small 



Vq + 4:Hvc 



(55) 
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indicates that such regime is incompatible with the boundary condition (13S|) . (14"3]) . 
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Figure 5. Growth of scalar perturbations in matter on the DGP brane is shown 
for different boundary conditions in the brane-bulk system. (The comoving 
wave number k/aoHo = 20 in all cases.) Boundary conditions are specified 
by (|38p and differ in the expression for A; namely: (a) the Koyama-Maartens 
condition A = —1/2, (b) A = —(1 + 37)/4, (c) the minimal condition ^4 = 0, 
and (d) A = 1/2. Condition (b) with time-dependent value of A was chosen 
because it simplifies equation (|27p. and condition (d) because it lies well inside 
the stability domain (|4ip . For comparison, note that A/a = 1 for all values of 
k in the standard CDM model with fi m = 1. 

The strong fc-dependence of the evolution of perturbations can be explained by the 
presence of the term V 2 t>c on the right-hand side of fl53|) . which leads to the generation 
of large perturbations of dark radiation 5pc- The quantity vq is being generated by the 
right-hand side of equation (1521) . The instability in the growth of perturbations for the 
Koyama-Maartens boundary condition is in agreement with the fact that the value of 
A = —1/2 lies well beyond the stability domain (H21 . 

As demonstrated earlier, depending upon the value of A, perturbations on the brane 
can be either unstable or quasi-stable. By unstable is meant A/a ^> 1 while quasi- 
stability implies A/a ~ 0(1). The quasi-stable region (14T!) is illustrated in figure HJ in 
which we show the results of a numerical integration of equations fl34l) - fl36l) for A = 1/2. 
It is instructive to compare this figure with the left panel of figure [2j One clearly sees the 
much weaker growth of perturbations as well as their scale-independence in this case. 
We therefore conclude that boundary conditions can strongly influence the evolution 
of perturbations on the brane. Our results are summarized in figure [5j which shows 
the evolution A/a obtained by integrating the system (13"l|) - (l3"t3|) for different boundary 
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conditions. (Results for the wave number k/aoH = 20 are shown.) We see that the 
growth of perturbations becomes weaker as the value of A approaches the stability 
domain (EH), and quasi-stability is observed for A = 1/2. 

The behaviour of scalar perturbations on the DGP brane in the case of scale-free 
boundary conditions (145j) is very similar to that shown in figure HI As in [HI [7] , we also 
find that perturbation growth on the DGP brane is slower than that in LCDM. 

5. Scalar perturbations in mimicry models 

5.1. Cosmic mimicry on the brane 

In our paper [IT] , we described a braneworld model in which cosmological evolution 
proceeds similarly to that of the Friedmannian cosmology but with different values of 
the effective matter parameter fl m , or, equivalently, with different values of the effective 
gravitational constant Gn, at different cosmological epochs. 

Specifically, for a spatially flat universe with zero background dark radiation 
(C = 0), the cosmological equation ( TlBl) can be written as follows: 



where p an d H are the energy density of matter and Hubble parameter, respectively, 
at the present moment of time. Proceeding from ({TBI) to this equation, we traded the 
brane tension a for the new parameters Hq and po- 

In the effective mimicry scenario [UJ, the parameters Hq—A/6 and l/£ 2 are assumed 
to be of the same order, and much larger than the present matter-density term po/3m 2 . 
The mimicry model has two regimes: one in the deep past (high matter density) and 
another in the more recent past (lower matter density). In the deep past, we have 




(56) 




(57) 



and the universe expands in a Friedmannian way 




(58) 



In the recent past, we have 




(59) 



and the expansion law is approximated by 




(60) 




(61) 
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is the parameter introduced in [TT]. In the case of the cosmological branch with upper 
sign, it is assumed that the coefficient a/ (a — 1) in fl60l) and in similar expressions is 
always positive, i.e., a is assumed to be greater than unity in this case. 

One can interpret the result (1601) either as a renormalization of the effective 
gravitational constant relative to its value in the deep past or as a renormalization 
of the effective density parameter: 

„LCDM _ „LCDM 

H 2 - Hi + P — £ , pLCDM = a 

6m z a =1= 1 

In a universe in which matter dominates in the energy density p, introducing the 

cosmological parameters 

3m 2 tf 2 ' l ~ PH%' 6tf 2 
one can express the Hubble parameter as a function of redshift for the case under 
consideration: 



fi m - „., rr2 , ^ - ^TTf ' fi A - _ Z772 ' ( 63 ) 



2/ 



2 



fi m (l + zr + 1 - tt m + 2fi £ =f 2a/^ \/l + 




±2J^WO m (l+z) 3 -O in + Jl + A TV^ • ( 64 ) 



For large values of the extra-dimensional parameters Q\ and Qe , for which the 
domain of redshifts z exists such that 

2 



f] m (i+ 2 ) 3 «^i + n AT ^j , (65) 

expression fl64|) in this domain reduces to 

^CDM (1 + z)3 + 1 _ n LCDM ; (66) 



H\z) 

tj2 '~ u "m V.- 1 - 1 ") 1 u "m 
tin 



l 



where 



a =F 1 v 
Equation fl66|) is equivalent to fl60|) . 

Remarkably, the behaviour of the Hubble parameter on the brane practically 
coincides with that in LCDM at low redshifts z <C z m , where 

, m = ^ (68) 

J tm 

This property was called 'cosmic mimicry' in [11] and z m the mimicry redshift. An 
important aspect of the mimicry model, illustrated by (16"B"j) and (16"7j) . is that the matter 
density entering fl66l) is an effective quantity. A consequence of this is the fact that 
the two density parameters Q m and f^ 001 ^ need not be equal and, based solely on 
observations of the coordinate distance, a low/high density braneworld could easily 
masquerade as LCDM with a moderate value f2^ CDM ~ 0.2 — 0.3. This situation is 
illustrated in figure El 
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Figure 6. This figure illustrates that -ffbranc — -f^LCDM for low redshifts, whereas 
at higher redshifts -ffbrane > -^LCDM (for Mimicry 1, shown in the left panel) 
or -ffbrane < ^LCDM (for Mimicry 2, shown in the right panel). This behaviour 
arises when fi m > fi^CDM (Mimicry 1) or ft m < fi^CDM (Mimicry 2). See 
[llj for details of mimicry models and for some cosmological consequences of 
mimicry. 



We remind the reader that the two signs in (156]) or (l64j) correspond to 
two complementary possibilities for embedding the brane in the higher-dimensional 
(Schwarzschild-AdS) bulk. In our ensuing discussion, we shall refer to the model with 
the lower (upper) sign in (1641) as Mimicry 1 (Mimicry 2), and we consider these two 
models separately. 

5.2. Cosmological perturbations in the mimicry model: 
I. Scale- dependent boundary conditions 

As expected, perturbations in the mimicry model crucially depend upon the type of 
boundary condition which has been imposed. Generally speaking, brane perturbations 
grow moderately for BC's which lie in the stability domain (14T1) or (1421) and more rapidly 
in the instability region. This remains true for mimicry models. In this section, we 
explore the behaviour of perturbations in this model for the boundary condition A = 0, 
which belongs to the instability class. In the next section, we shall explore BC's which 
give rise to more moderate and scale-independent behavior. 

The growth of perturbations if A = is substituted in (1551) is illustrated in figure CD 
The fc-dependence, clearly seen in this figure, can be understood by inspecting the 
system of equations (I34|) - (l36|) . Even if we start with zero initial conditions for the dark- 
radiation components Spc and vc, the nontrivial right-hand side of Eq. ( 135|) leads to the 
generation of vq\ then, via the /^-dependent right-hand side of (|36|) . the density Spc is 
generated, which later influences the growth of perturbations of matter in ( }34"1) . The 
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instability in the growth of perturbations is explained by the fact that the value of A = 
lies outside the stability domain (|4ip or (T421) . However, the growth of perturbations is 
not as dramatic in this case as in the DGP model with the Koyama-Maartens BC's, 
mainly because the value A = lies much closer to the boundary ([4*2]) than the Koyama- 
Maartens value A = —1/2. 
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Figure 7. Growth of matter perturbations in the Mimicry 1 and Mimicry 2 
models are shown for different values of the comoving wave number h/a^Ho 
(indicated by numbers above the corresponding curves) and for the minimal 
boundary condition A = in (|38p . Both models have the same effective 
parameter fi^CDM = 0.22, hence, quite different matter content, indicated by 
the parameter f2 m . The value of vq is set to zero initially. The position of the 
mimicry redshift z m is indicated by the vertical dotted lines. 



Qualitatively, the evolution of matter perturbations in mimicry models can be 
understood as follows: during the early stages of matter-domination the last term on 
the right-hand side of equation (1511) is not very important, which transforms (1511) into 
a closed equation for the matter perturbation. Indeed, in the pre-mimicry regime, for 
z ^> z m , we have \i 2 X\ ^> 1 for the quantity in the denominator of the last term 
on the right-hand side of ( 1341) . which makes this term relatively small for moderate 
values of 5pc- Thus, perturbations in matter evolve according to (|46l) on all spatial 
scales, for redshifts greater than the mimicry redshift z m . For z < z m , the quantity 
|£ 2 A| is of order unity. By this time, the perturbations Spc have grown large, and their 
amplitude strongly depends on the wave number. Through the last term in equation 
( I34p . they begin to influence the growth of matter perturbations for z ~ z m , resulting 
in the /c-dependent growth of the latter. The reason for the opposite fc-dependence of 
matter perturbations in Mimicry 1 and Mimicry 2 shown in figure [7] is connected with 
the difference in the sign of A — defined in (132]) — for the two models. Thus, the last 
term in (1341) comes with opposite signs in Mimicry 1 and Mimicry 2, and therefore works 
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in opposite directions in these two models. 

Well inside the mimicry regime, for z <C z m , we have 7 ~ 1/3, so that the second 
term on the right-hand side of ( 1351) can be ignored if matter perturbations are not too 
large. Then equations (|35|) . (1361) . and (I38I) lead to a closed system of equations for the 
evolution of dark-radiation perturbations. Substituting 7 = 1/3 into this system, we 
obtain: 



where we assumed A 56 —1/2 to be constant. The function £ obeys an oscillator-type 
equation 



This means that both Spc and vq rapidly decay during the mimicry regime (oscillating 
approximately in opposite phase) and the last term on the right-hand side of (134]) again 
becomes unimportant. In particular, this will describe the behaviour of the mimicry 
model with the minimal boundary condition A = 0. The transient oscillatory character 
of 5pc induces transient oscillations with small amplitude in A through the last term in 
( 136]) . These small oscillations can be noticed in figure [7] for log (1 + z) > 1, particularly 
for values 40 and 30 of the comoving wave number |+] 

Two important features of mimicry models deserve to be highlighted: 

(i) As demonstrated in figure [7J there is a strong suppression of long- wavelength modes 
in Mimicry 2. 

(ii) From this figure, we also find that the growth of short-wavelength modes in 
Mimicry 2 can be substantial, even in a low-density universe. 

Both properties could lead to interesting cosmological consequences. For instance, 
the relative suppression of low- A; modes may lead to a corresponding suppression of 
low-multipole fluctuations in the cosmic microwave background, while the increased 
amplitude of high-fc modes could lead to an earlier epoch of structure formation. (Since 
the mimicry models behave as LCDM at low redshifts, they satisfy the supernova 
constraints quite well.) A detailed investigation of both effects, however, requires that 
we know the form of the transfer function of fluctuations in matter (and dark radiation) 
at the end of the radiative epoch. Such an investigation lies outside the scope of the 
present paper, but we may return to it elsewhereEl 

+ For the Koyama-Maartens boundary condition A = —1/2, the approximation described above is not 
valid during the mimicry stage. Instead, during mimicry, the value of vq decays without oscillating 
approximately as vq oc 1/a 3 , as can be seen from equations (|35[) . (|38[) . and the value of Spc also decays, 
which follows from ([36)) . 

* For simplicity, the amplitudes of all fc-modes were assumed to be equal at high redshifts in figures 
[31 [5] and [7] A more realistic portrayal of A(fc) should take into consideration the initial spectrum and 
the properties of the transfer function for matter and dark radiation, and we shall return to this in a 
future work. 
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For the minimal boundary condition (A = 0), assumed in this section, equations 
(]28j) and (PJ imply $ = ^ and 



which is a generalization of the Poisson equation for the mimicry brane. 

Mimicry models with the Koyama-Maartens boundary condition exhibit much 
stronger instability in the growth of Spc for high values of k, enhancing the growth 
of matter perturbations (not shown) . This can be explained by the fact that A = is 
much closer to the boundary of the stability domain fH2l) than the Koyama-Maartens 
value A = —1/2. For the latter, the density perturbation A in the Mimicry 1 model 
grows to be large and negative, while the perturbation Spc becomes large and positive; 
for instance, both A and the dimensionless quantity 5pc/m 2 H 2 grow by a factor of 
10 11 for k/aoHo = 40. This provides another example of the very strong dependence of 
perturbation evolution on boundary conditions. 

In our calculations, we have not found any significant dependence of the eventual 
growth of perturbations on initial conditions for dark radiation specified in a reasonable 
range (at z = 10 3 ). 

5. 3. Cosmological perturbations in the Mimicry model: 
II. Scale-free boundary conditions 

As mentioned earlier, BC's lying in the stability region ( T4T|) lead to an almost scale-free 
growth of density perturbations. A similar result is obtained if we assume the scale-free 
boundary condition ()45l) of section 13.21 In this case, the momentum potential vc decays 
as Vc oc a -4 , and its spatial gradients in fl36|) can therefore be neglected. The same is 
true, of course, if one considers super-horizon modes with k <C aH. In both cases, we 
have approximately 5pc oc 1/a 4 , suggesting that the dynamical role of perturbations in 
dark radiation is unimportant. This results in a radical simplification: as in the DGP 
model, for BC's lying in the stability region, the growth of perturbations in matter can 
be effectively described by a single second-order differential equation (T4l)]) . namely, 



where A and 7 are defined in (13~TT) and (|32|) . respectively. We shall call Q(z) in (!72j) 
the 'gravity term' since it incorporates the effects of modified gravity on the growth of 
perturbations. The value of this term on the brane can depart from the canonical = 1 
in general relativity. 

Figure M shows the behaviour of Q(z) for a typical mimicry model. At redshifts 
significantly larger than the mimicry redshift, z ^> z m , we have Q(z) ~ 1, whereas at 
low redshifts, z <C z m , the value of Q(z) changes to 




(71) 



A + 2HA = 6 





(72) 



,LCDM 



e(z) 



p 



for z z, 
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0(z) 



Mimicry on Brane 1 



1 2 
Loq( 1 +z) 



Figure 8. The Hubble parameter in the braneworld 'Mimicry 1' is shown 
relative to LCDM (solid). The LCDM model has ^cdm = 0.22 while 
= 0.27 on the brane. Also shown is the 'gravity term' Q(z) defined in 
(|46p whose value diminishes from unity at high redshifts to the asymptotic 
form (|73p at low redshifts. The dotted vertical line shows the mimicry redshift 

Zm. ~ 37. 



where 

p LCDM ig 

defined in f[6"2l) . The solid line in the same figure shows the ratio 
of the Hubble parameter on the brane to that in LCDM. The consequences of this 
behaviour for the growth equation (1151) are very interesting. Substituting ([75]) into (14"61) 
and noting that H(z) ~ H LCDM during mimicry, we recover the standard equation 
describing perturbation growth in the LCDM model 

JjCDMa 

(74) 



A + 2H 



LCDM 



A 



2m 2 



Thus, ordinary matter in mimicry models gravitates in agreement with the effective 
value of the gravitational constant which appears in the cosmological equations (1581) 
and (ISP]) . 

We therefore conclude that, deep in the mimicry regime (z <C z m ), perturbations 
grow at the same rate on the brane and in LCDM. This is borne out by figure [9] 
which shows the results of a numerical integration of ( 146]) for Mimicry 1 [integrating the 
exact system (I34l) - (l36j) gives indistinguishable results]. Notice that the total amplitude 
of fluctuations during mimicry in this model is greater on the brane than in LCDM. 
Indeed, for mimicry models, we have 

Abrane 



A 



LCDM 



J^LCDM 



for z -C z n 



(75) 
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Figure 9. Density perturbations on the Mimicry 1 brane (dot-dash) and in 
LCDM (solid). The evolution of perturbations in Mimicry 1 in this case is 
effectively described by (|46p . In both cases, the perturbation amplitude is scaled 
by the expansion factor a(t). (It may be noted that A/a = 1 in standard CDM 
with O m = 1.) The dotted vertical line shows the mimicry redshift z m ~ 37. 
The braneworld has J7 m = 0.27 while f2 m CDM = 0.22. This leads to a moderate 
enhancement in the amplitude of brane perturbations over LCDM. 



and this ratio is greater than unity for Mimicry 1 . 

Since the contribution from perturbations in dark radiation can be neglected, the 
growth of matter perturbations in Mimicry 2 is again described by (Hoi) and by f!75|) . 
However, since fl m < f2 m CDM in this case, the final amplitude of perturbations will be 
smaller in Mimicry 2 than the corresponding quantity in LCDM, which is the opposite 
of what we have for Mimicry 1 . 

It is interesting that during mimicry, when 7 w 1/3, the relation between the 
gravitational potentials $ and \I/ reduces to the general- relativistic form $ = as can 
be seen from (jlS|) . where $ satisfies the generalized Poisson equation (jl7|h namely, 

lv 2 $ = 6^L . (76 ) 

An interesting feature of Mimicry 1 is that, at early times, the expansion rate in this 
model exceeds that in LCDM, i.e., H (z) iMimiayi > -^(^)|lcdm for z > z m (see Figs. |6]& 
|8]). [The opposite is the case for Mimicry 2: the expansion rate in this model is lower 
than that in LCDM at early times, i.e., H(z)\ Mimicry 2 < -^(-z)Ilcdm for z > z m .} As we 
can see, this has important consequences for the growth of structure in this model. The 
increase in the growth of perturbations in Mimicry 1 relative to LCDM occurs during 
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the period before and slightly after the mimicry redshift has been reached, when the 



not reached its asymptotic form ( 173]) . A lower value of H(z) in (06]) diminishes the 
damping of perturbations due to cosmological expansion while a slower drop in Q(z) 
signifies a much more gradual decrease in the force of gravity. Consequently, there is a 
net increase in the growth of perturbations on the Mimicry 1 brane relative to LCDmJ^I 
For the models in figure [9] which have f2j^ CDM = 0.22 and Q m = 0.27, the increase is 
about 20%. The increased amplitude of perturbations in Mimicry 1 stands in contrast to 
the DGP model as well as Quintessence model, in both of which linearized perturbations 
grow at a slower rate than in the LCDM cosmology P, [7J, [29] . 

It is important to note that observations of galaxy clustering by the 2dFGRS survey 
provide the following estimate [30] for perturbation growth at a redshift z = 0.15 : 



where 5 = 5p/p. Since the growth of perturbations during the mimicry regime stays 
proportional to that in the LCDM model (<5 mimic oc <5lcdm, z ^ z m), it follows that if 
perturbations in the LCDM model satisfy ( 1771) (which they do), then so will those in 
the mimicry scenario. Nevertheless, as we have seen, the net increase in the amplitude 
of perturbations on the brane is larger than that in the LCDM model. This clearly has 
important cosmological consequences since it can enhance structure formation at high 
redshifts as well as lead to higher values of as- Thus, while preserving the many virtues 
of the LCDM model, the mimicry models add important new features which could be 
tested by current and future observations. We hope to return to some of these issues in 
a future work. 

6. Discussion 

Braneworld theories with large extra dimensions, while having a number of very 
attractive properties, also have a common difficulty: on the one hand, the dynamics of 
the higher- dimensional bulk space needs to be taken into account in order to understand 
brane dynamics; on the other hand, all observables are restricted to the four- dimensional 
brane. In field-theoretic language, the situation can be described in terms of an infinite 
(quasi) -continuum of Kaluza-Klein gravitational modes existing on the brane from 
the brane viewpoint. This property makes braneworld theory complicated, solutions 
on the brane non-unique and evolution nonlocal. For instance, while the solution 
with a spherically symmetric source is relatively simple and straightforward in general 
relativity, a similar problem in braneworld theory does not appear to have a unique 
solution (see, e.g., [3~T]). 

jj Figure [5] clearly shows that H(z) reaches its asymptotic form much sooner than Q(z). Notice that, 
at redshifts slightly larger than z m , the value of Q(z) exceeds unity. The dependence of perturbation 
growth on the mimicry redshift z m is very weak, and (|75p is a robust result which holds to an accuracy 
of better than 2% for a wide range of parameter values. 



relative expansion rate H(z)/H " LCDM is declining while the 'gravity term' Q(z) has still 




0.51 ±0.11, 



(77) 
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Fortunately, in situations possessing a high degree of symmetry, the above 
properties of braneworld theory do not affect its cosmological solutions (at least, in 
the simplest case of one extra dimension). Thus, homogeneous and isotropic cosmology 
on the brane is almost uniquely specified since it involves only one additional integration 
constant which is associated with the mass of a black hole in the five- dimensional bulk 
space. This makes braneworld theory an interesting alternative for modelling dark- 
energy [22] and dark-matter [11] effects on cosmological scales. 

However, in order to turn a braneworld model into a complete theory of gravity 
viable in all physical circumstances, it is necessary to address the issue of boundary 
conditions. Usually, one tries to formulate reasonable conditions in the bulk by 
demanding that the bulk metric be nonsingular or by employing other (regulatory) 
branes. However, neither of these conditions have been implemented in braneworld 
theory in full generality; moreover, they leave open the problem of nonlocality of the 
gravitational laws on the brane (since the brane is left open to influences from the bulk). 

In this paper, we adopted a different approach to the issue of boundary conditions 
in the brane-bulk system. From a broader perspective, boundary conditions can be 
regarded as any conditions restricting the space of solutions. Our approach is to specify 
such conditions directly on the brane which represents the observable world, in order to 
arrive at a local and closed system of equations on the brane. The behaviour of the metric 
in the bulk is of no further concern in this approach, since this metric is for all practical 
purposes unobservable. Since the nonlocality of the braneworld equations is known to 
be connected with the dynamical properties of the bulk Weyl tensor projected onto the 
brane [15], it is natural to consider the possibility of imposing certain restrictions on 
this tensor. Perhaps, the simplest condition is to set to zero its (appropriately defined) 
anisotropic stress. This is fully compatible with all the equations of the theory and 
results in a brane universe described by a modified theory of gravity and having an 
additional invisible component — dark radiation — which is endowed with nontrivial 
dynamics. More generally, we suggest (ITTj) and (|38|) which describes a one-parameter 
family of BC's with the parameter A. This family generalizes the boundary condition 
derived by Koyama and Maartens [7] for the DGP model in the small-scale and quasi- 
static approximation (Koyama and Maartens derived the value A = —1/2). 

An important conclusion of our paper is that the growth of perturbations in 
braneworld models strongly depends upon our choice of BC's. This was illustrated 
in figure [5] for the DGP model. Specifying boundary conditions in the form ( 1381) allows 
us to determine regions of stability and instability in terms of the single parameter 
A; they are described by Eq. f j4TI) . In the DGP model, perturbations are explicitly 
demonstrated to be quasi-stable for A = 1/2 (figure H]) and unstable for A = —1/2 
(figure [2]). In the instability domain, gradients in the momentum potential Vq of dark 
radiation, lead to the creation of perturbations in this quantity via equation (1351) . This 
effect can significantly boost the growth of perturbations in matter. An important 
implication of this effect is that perturbations in the baryonic component might overcome 
the 'growth problem,' which plagues them in standard general relativity, and grow 
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to acceptable values without requiring the presence of (deep potential wells in) dark 
matter. The Mimicry 2 model looks promising from this perspective. Note that, in this 
model, the expansion of a low-density universe is virtually indistinguishable from that 
of a (higher-density) LCDM model. Enhanced perturbation growth in this 'mimicry' 
model might permit a low-density universe to also match observations of gravitational 
clustering (thereby circumventing the need for large amounts of dark matter to achieve 
this purpose) Jffi A detailed investigation of this effect, however, lies outside the scope 
of the present paper and will be taken up elsewhere. 

Values of A lying in the stability region fHTj) or a scale-free boundary condition 
such as (H5l) may also be important. In this case, perturbations in the DGP model grow 
slower than in LCDM whereas perturbations in Mimicry 1 grow somewhat faster. This 
suggests that structure formation may occur slightly earlier in Mimicry 1 than it does 
in LCDM. 

Two related points deserve mention. In general relativity, the linearized 
perturbation equation 

5 + 2H5 = lH^n m (^-) 3 5 (78) 



2 u V a 

depends exclusively upon the expansion history of the universe and the value of Q m , so 
that a knowledge of H(z) from independent sources (such as the luminosity distance) 
permits one to determine 5(z) and, vice-versa, the observed value of 5(z) can be used 
to reconstruct the expansion history H(z) [51] . As we have demonstrated in this paper, 
equation (1781 is no longer valid for the braneworld, and it remains to be seen whether 
knowing H(z) and Q m will allow us to reconstruct 5(z) uniquely. 

It is also well known that the expansion history, H(z), does not characterize a 
given world model uniquely, and it is conceivable that cosmological models having 
fundamentally different theoretical underpinnings (such as different forms of the matter 
Lagrangian or different field equations for gravity) could have identical expansion 
histories (some examples may be found in [31] which also contains references to earlier 
work). As an illustration, consider the Brane 2 model in which the effective equation 
of state is never phantom-like, so that w e f[ > — 1. As we have seen, Mimicry 2 and the 
DGP model form important subclasses of this braneworld. A quintessence potential 
corresponding to a given expansion history can be determined from [35] 

8tt<7 . . H 2 x dH 2 1 3 

*v<?) =^-^T-2^' (79) 

2 dlnH QmX x = l + z. (80) 




3H^x dx H 2 

Substitution for H(z) from ({TBI) with the plus sign results in a quintessence model which 
has the same expansion history as the braneworld. Perturbative effects, however, are 

ff Of course, a valid model of structure formation should also agree with observations of gravitational 
lensing which appear to require some amount of dark matter 32J , as well as provide a good fit to the 
observed power spectrum of galaxies and fluctuations in the cosmic microwave background; see [33] for 
discussions of this issue within a different context. 
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likely to break this degeneracy, since it is unlikely that matter perturbations will grow 
at the same rate in the quintessence model as on the brane. 

Although we have explicitly solved the perturbation equations only for two cases: 
the mimicry models [11] and the DGP braneworld [2], it is quite clear that the treatment 
developed by us may have larger ramifications since it could be applied to a broader 
class of models. In [36], for instance, it was shown that a braneworld dual to the 
Randall-Sundrum II model - - the bulk dimension being time-like instead of space- 
like - - will bounce at early times, thereby avoiding the singularity problem which 
plagues general relativity. It would be an interesting exercise to apply the formalism 
developed in this paper to this non-singular braneworld to see if the problem of tachyonic 
gravitational modes can be avoided. Since the boundary conditions of type (fTU]) - ffTT|) . 
specified directly on the brane, effectively freeze certain degrees of freedom in the brane- 
bulk system, the corresponding braneworld theory may also be free from the debatable 
problem of ghosts (see [37] for a recent discussion of this problem). 

Whether boundary conditions such as those described by equations ( fT0l) - (|T7l) or 
f|45|) will remain in place for a more fundamental extra-dimensional theory is presently a 
moot point. Perhaps, by comparing the consequences of different boundary conditions 
with observations we will gain a better understanding of the type of braneworld theory 
most consistent with reality. 

We believe that our calculations have opened up an interesting avenue for further 
research in braneworld cosmology. Related astrophysical phenomena which demand 
further exploration include gravitational lensing, fluctuations in the cosmic microwave 
background, the power spectrum of density fluctuations, etc. For a deeper understanding 
of gravitational instability on the brane, we also need to examine the behaviour of 
perturbations during the radiative epoch as well as their origin as vacuum fluctuations 
(presumably) during inflation. Additional issues of considerable interest relate to non- 
linear gravitational clustering on the brane. For instance, it is well known that a density 
field which was originally Gaussian develops non-Gaussian features even in the weakly 
non- linear regime of gravitational instability [3H]. Although these effects do not appear 
to be unduly sensitive to the presence of dark energy, this has only been tested for 
ACDM and a few other conventional dark-energy models [39] . The value of the lower 
moments of the density field such as its skewness (<5 3 ) / (5 2 ) 2 and kurtosis (5 4 ) / (5 2 ) 3 in 
models which give rise to acceleration because of amendments to the gravitational sector 
of the theory remain unexplored and could lead to a deeper understanding of structure 
formation in these models. 

We hope to return to some of these issues in the future. 
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Appendix A. Attempts at deriving boundary conditions 

Recent attempts at deriving proper boundary conditions for the brane-bulk system by 
considering perturbations in the bulk were made in [3 [8] . Here we would like to discuss 
this approach in more detail and illustrate some difficulties with this treatment. 

Both approaches study perturbation in the DGP model, which has Minkowski bulk 
background, and use the Mukohyama master equation [9] for the scalar variable Q(t, x, y) 
written in Gaussian normal coordinates with respect to the brane, which is located at 
y — in the fifth coordinate y and has intrinsic coordinates {t, x}: 

di \nb* dt ) + dy [b 3 dy ) ' ( } 

where 

n(t, y) = l+(^ + H ) jy, bit, y) = a(l + Hy) (A.2) 

are the functions describing the unperturbed bulk metric 

ds 2 = -n 2 dt 2 + b 2 dx 2 + dy 2 , (A.3) 

and k is the wave number in the Fourier decomposition in the x space. 

The (suitably normalized) Mukohyama variable is related to perturbations on the 
DGP brane via ^U\: 



i 



6pc = ~h n{ y=° ' 




v c = —s [ Q ~ HQ ) y=0 > (A.4) 



y=o 



i2 f)0 s 

Stt c = -— ( 3Q - 3HQ + — Q + 3H (l - E~ x \ — 
6a \ a 2 v oy 

Koyama and Maartens [7j study perturbations on subhorizon scales, k 3> aH, and 
employ what they call 'quasi-static approximation,' which consists in neglecting the time 
derivatives in flA.ll) and flA.4j) . By considering the approximate 'quasi-static' solution 
for Q with appropriate (decaying) boundary condition in the bulk, they argue that the 
term with the derivative dVt/dy in the brackets on the right-hand side of the last relation 
in ( 1A.4I) can also be neglected, and this relation then becomes 

6n c ~ Q\ y=Q . (A.5) 

This results in the Koyama-Maartens boundary condition 
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However, as was shown in the paper, the boundary condition (1551) . (l4"3"j) leads 
directly to equation (HOl) . which exhibits strong scale-dependent instability The last 
effect was also demonstrated numerically. We note that, in deriving equation ( HUI) . we 
used the exact system of equations ([34")) - ([36]) on the brane, and we did not perform 
any time differentiation of the approximate relation (1381) . which might be illegitimate 
in the quasi-static approximation. This means, in our view, that the quasi-static 
approximation, which was assumed in deriving conditions f)A.5j) . cannot hold during 
the whole course of evolution described by the system of equations for perturbations 
(E1-(J36D. 

Sawicki, Song, and Hu [8] study scalar perturbations on all scales. They note that 
the distance y^ to the causal horizon of the brane remains proportional to H~ x during 
periods of domination of matter with constant equation of state, and they impose the 
zero boundary condition for Q at this causal horizon. This does not fix the value of Q 
at the brane, and it is done separately by setting a 'scaling' ansatz of the form 

Q = A(p)a p G (x, 4t J , * = — , 2/hor = £H~ X . (A.6) 

They substitute (1 A. 61) into ( 1A.1I) to derive a differential equation on G (x, during 
periods of domination of a fluid with constant equation of state (radiation or matter). 
In doing this, they treat p, A(p), and £ as constants and also neglect the derivative 
of G (x, with respect to the second argument. One cane note that this already 
imposes certain restriction on the form of Q on the brane, hence, via (1A.4I) . also on the 
perturbations of dark radiation. The result is the differential equation for G (x, -^j of 
the form 

fo{p,x) - I \- ) g{p,x) 



d 2 G | dG 
dx 2 ' dx 



G = 0, (A.7) 



where fo, fx, and g are some functions, and in which the quantity -kr plays the role of a 
parameter. This equation is solved with the boundary conditions G(0) = 1, G(l) = 0, 
which uniquely determines the derivative G'(0), hence, the quantity dfl/dy on the brane 
needed in ( 1A.4I) . 

The power p in (1A.61) is not yet specified, and it is adjusted in [S] by an iterative 
procedure so as to satisfy the Bianchi identity on the brane coupled to the cosmological 
equations for matter perturbations. In this way, p also becomes a function of the scale 
factor a. 

One can see that a number of assumptions were made in this procedure. In deriving 
equation (IA.7I) . the time derivatives of p and the partial derivative of G (x, -^j with 
respect to its second argument were neglected. Together with ansatz (IA.6I) . this means 
that an effective 'quasi-static' approximation was used for Q, similar to the one employed 
in [7]. It is not surprising, therefore, that, on subhorizon scales, the result of [7] is 
reproduced, with the approximate relations (13"%|) . (14"3"|) . for which the quasi-static regime 
is unstable, as we have shown in this paper. 
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Furthermore, the authors did not extend smoothly their approximate solution in 
the bulk beyond the horizon y = y\ lov , so it is not at all clear whether it is really 
singularity-free in the bulk, as was claimed. To specify simply Q = for y > y^ or 
means introducing a discontinuity in the derivative of at y — ytor, which itself can 
be regarded as a singularity. This once again illustrates the general difficulty which one 
encounters in finding singularity-free solutions in the bulk. 
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